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^ The recent observations support that our Universe is flat and expanding with acceler- 

ation. We analyze a general class of quintessence models by using the recent type la 
supernova and the first year Wilkinson Microwave Anisotropy Probe (WMAP) obser- 
vations. For a flat universe dominated by a dark energy with constant w which is a 
special case of the general model, we find that f2 m o = 0.30^gQg and uiq < —0.82, and 
the turnaround redshift zt when the universe switched from the deceleration phase to 
the acceleration phase is z-p = 0.65. For the general model, we find that f2 m o ~ 0.3, 
W Q0 ~ —1-0, /3 ~ 0.5 and zt ~ 0.67. A model independent polynomial parameterization 
^ of dark energy is also considered, the best fit model gives f2 m o = 0.40 ±0.14, loqq = — 1.4 

. and zt = 0.37. 
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a: 

1. Introduction 

The type la supernova (SN la) observations indicate that the expansion of the 
Universe is speeding up rather than slowing down - . The measurement of the 
5— i ' anisotropy of the cosmic microwave background (CMB) favors a fiat universe 5 > 6 ' 7 . 

The observation of type la supernova SN 1997ff at z ~ 1.7 also provides the evidence 
that the Universe is in the acceleration phase and was in the deceleration phase in 
the past 4,8 . The transition from the deceleration phase to the acceleration phase 
happened around the redshift zt ~ 0.4 4:9 . In this paper, we use the notation zt for 
the transition redshift. A new component with negative pressure widely referred as 
dark energy is usually introduced to explain the accelerating expansion. The sim- 
plest form of dark energy is the cosmological constant with the equation of state 
parameter u\ — — 1 . One easily generalizes the cosmological constant model to dy- 
namical cosmological constant models such as the dark energy model with negative 
constant equation of state parameter — 1 < luq < — 1/3 and the holographic dark 
energy models 10 . If we remove the null energy condition restriction luq > — 1 to 
allow supernegative luq < — 1, then we have the phantom energy models 11 . More 
exotic equation of state is also possible, such as the Chaplygin gas model with the 
equation of state p = —A/p and the generalized Chaplygin gas model with the 
equation of state p = —A/p a 12 . In general, a scalar field Q that slowly evolves 
down its potential V(Q) takes the role of a dynamical cosmological constant. The 
scalar field Q is also called the quintessence field 13 - 20 . The energy density of the 
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quintessence field must remain very small compared with that of radiation or mat- 
ter at early epoches and evolves in a way that it started to dominate the universe 
around the redshift 0.4. Instead of the quintessence field with the usual kinetic term 
Q 2 /2, tachyon field as dark energy was also proposed 21 . The tachyon models have 
the accelerated phase followed by the decelerated phase. 

Although most dark energy models are consistent with current observations, the 
nature of dark energy is still mysterious. Therefore it is also possible that the obser- 
vations show a sign of the breakdown of the standard cosmology. Some alternative 
models to dark energy models were proposed along this line of reasoning. These 
models arc motivated by extra dimensions. In these models, the usual Friedmann 
equation H 2 — 8irGp/3 is modified to a general form H 2 = g(p) and the universe 
is composed of the ordinary matter only 22 - 29 . In other words, the dark energy 
component is unnecessary. 

In this paper, I first use the 58 SN la data in Ref. 3, the 186 SN la data in 
Ref. 4 and WMAP data 7 to constrain the parameter space of a general class of 
quintessence models discussed in Ref. 20. In that model, a general relation between 
the potential energy and the kinetic energy of the quintessence field was proposed. 
As we know, the average kinetic energy is the same as the average potential energy 
for a point mass in a harmonic oscillator. For a stable, self-gravitating, spherical 
distribution of equal mass objects, the total kinetic energy of the objects is equal 
to minus 1/2 times the total gravitational potential energy. Therefore, the physics 
of dark energy may be determined if the relationship between the potential energy 
and the kinetic energy is known. Then I consider three different model independent 
parameterizations of loq to find out some properties of dark energy. After we de- 
termine the parameters in these models, the transition redshift zt is obtained. The 
paper is organized as follows. After a brief introduction in section 1, the general 
class of models is reviewed in section 2. In section 3, 1 discuss the methodology used 
in this paper. In section 4, 1 give the main fitting results. In section 5, 1 conclude the 
paper by using a model independent analysis and compare the results with those 
in the literature. 



2. Model Review 

For a spatially flat, isotropic and homogeneous universe with both an ordinary pres- 
sureless dust matter and a minimally coupled scalar field Q source, the Friedmann 
equations are 

^g) 2 = ^(p m + P Q ), (1) 
a AttG , n . 

- = ^-(Pm + PQ +3p Q ), (2) 

Q + 3HQ + V'(Q) = 1 (3) 
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where dot means derivative with respect to time, p m = p m o(ao/ a ) 3 is the matter 
energy density, a subscript means the value of the variable at present time, pq — 
Q 2 /2 + V(Q), pq = Q 2 /2 - V(Q), V'(Q) = dV{Q)/dQ and V(Q) is the potential 
of the quintessence field. In Ref. 20, a general relationship 

V(Q) = 0Q 2 + C, (4) 

was proposed instead of assuming a particular potential for the quintessence field 
or a particular form of the scale factor, where and C are constants. Note that 
the above equation (4) is a constraint equation, one should not just substitute the 
above equation into the Lagrangian and thinks that the model is equivalent to 
a 1/2 + kinetic term plus a cosmological constant term C. The above general 
potential includes the hyperbolic potential and the double exponential potential. 
In terms of pqq and wqo, we have 
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To make the quintessence field sub-dominated during early times, we require that 
> 0.5. The transition from deceleration to acceleration happens when the decel- 
eration parameter q = —dH 2 /a = 0. From equations (2), (6) and (7), in terms of 
the redshift parameter 1 + z = a a /a, we have 

(1 + z T ) 3 + 2(1 - 0) (1 + c Q0 ) ^ (1 + z T ) 6/(2/3+1) - [1 - 20 - (1 + 20 Vqq] ^£ = 0.(9) 

This equation gives a relationship between wq and Oqo. Now let us turn our 
attention to two special cases. 

Case 1: C = 0, the equation of state of the scalar field is a constant, cjq = 
(1/2 - 0)/(l/2 + 0). The potential is 17 < 18 

V(Q) = A[smhk(Q/a + B)]- a , 

where a = 2/(0 - 1/2), k 2 = 48irG/(20 + 1), A^ 1 ' 2 = (1/2 + 
0)Cf +1 00- 1 / 2 /(p mO a 3 o ) and B is an arbitrary integration constant. 

Case 2: — 1/2, the pressure of the scalar field becomes a constant pq = — C = 
wqoPqo and the potential is the double exponential potential 19 



A 2 

V(Q) = — [cxp(2aQ) + exp(-2aQ)] + 



6ttG 



A^_ 

4 ' 



The constant pressure model is equivalent to an ordinary matter with effective 
matter content f2£f — (1 + u)q )Qqo + ^mo phis a cosmological constant p\ = 
-^QoPqo- 
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3. Methodology 

In order to use the WMAP result, one usually parameterizes the location of the 
m-th peak of CMB power spectrum as 30 

l m = (TO - 4>m)lA, 

where the acoustic scale I a is 

7T T - T ls _ 7T J Q ZlS dzj yfgjz) 



Ia 



C s Tig 



the conformal time at the last scattering n s and at today t are 
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Jo Jo 
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aoH \/g(z)' 

g(z) = n m0 (l + z) 3 + ft r0 (l + ^) 4 + (1/2 + /3)(1 + ^ QO )r!Qo(l + zf^ 2fj+1 '> 
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+ /3 WQO 
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(10) 

(11) 
(12) 

(13) 



f2 r o = 8.35 x 10 5 is the current radiation component and z\ s — 1089 ± 1 6 . The 
difficulty of this method is that there are several undetermined parameters, such 

1 /2 

as (f) m and c s . Instead, we use the CMB shift parameter 1Z = a H (To — t\ s ) = 
1.716 ± 0.062 31 to constrain the model. 
The luminosity distance is defined as 



a (\ n ^ f dt ' 1 + z f 
d L (z) = o c(l + z) / —— = c— — 

Jt a{t ) H Jo 
The apparent magnitude rcdshift relation becomes 



du[g(u)} 



-1/2 



(14) 



m(z) = M + 51og 10 
= M + 5 log 10 



d L (z) + 25 = Ai + 51og 10 V L (z) 

c(l + z) f dulgiu)]- 1 ' 2 
Jo 



(15) 



where T) L (z) = H d L (z) is the " Hubble-constant-free" luminosity distance, M is 
the absolute peak magnitude and M. = M — 5 log 10 Ho + 25. M. can be determined 
from the low redshift limit at where T>l(z) — z. We use the 54 SNe Ia data with 
both the stretch correction and the host-galaxy extinction correction, i.e., the fit 
3 supernova data in Ref. 3 (we refer the data as Knop sample), and the 186 SNe 
Ia data in Ref. 4 (we refer the data as Riess sample) to constrain the model. The 
parameters in the model are determined using a ^-minimization procedure based 
on MINUIT code. There are four parameters in the fit: the current mass density 
f2 m o, the current dark energy equation of state parameter wq , the constant /3 as 
well as the nuisance parameter M. The range of parameter space is fi m0 = [0, 1] 
and cjqo = (-1, 0]. 
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4. Results 

For the dark energy model with constant u)q, i.e., the model with C = 0, the best 
fit parameters to the 54 knop sample are fi m o = [0, 0.46] centered at 0.30 and 
uiQ = (—1.0, —0.44] centered at almost —1.0 with x 2 = 45.6 at 68% confidence 
level. The best fit parameters to the 157 Riess gold sample are f2 m0 = [0.19, 0.37] 
centered at 0.31 and uiq = (—1, —0.75] centered at almost —1.0 with x 2 = 177.1 at 
68% confidence level. The best fit parameters to the 186 Ricss gold and silver sample 
are fi m o = [0.22, 0.36] centered at 0.30 and wq = (—1, —0.81] centered at almost 
— 1.0 with x 2 = 232.3 at 68% confidence level. The best fit parameters to the Riess 
gold sample and WMAP data combined are O m o = [0.22, 0.36] centered at 0.30 and 
ujq = (—1.0, —0.81] centered at almost —1.0 with x 2 = 177.3 at 68% confidence 
level. The best fit parameters to the Riess gold and silver sample and WMAP data 
combined are fi m0 = [0.23, 0.35] centered at 0.30 and ujq = (—1.0, —0.84] centered 
at almost —1.0 with x 2 — 232.4 at 68% confidence level. From the above results, it 
is shown that the addition of Riess silver data give almost the same results as those 
from Riess gold sample only. Therefore, we will use Riess gold sample only in the 
following discussion. The confidence regions of Sl m o and loq are shown in figure 1. 

For the general model C ^ 0, we first analyze the special constant pressure 
model (3 — 1/2 which is equivalent to the A-CDM mdoel. The best fits to the 54 
Knop sample are f2 m o = [0, 0.46] centered at almost zero and wqo = (—1, —0.54] 
centered at —0.71 with x 2 — 45.6 at 68% confidence level. The best fits to the 157 
Riess gold sample are: < f2 m o < 0.37 and — 1 < ojq < —0.63 at 68% confidence 
level with x 2 = 177.1. Note that the effective fi m0 ~ 0.3 although the best fit U m0 is 
almost zero. The best fits to the 157 Riess gold sample and WMAP data combined 
are: 0.25 < fi m0 < 0.35 and -1 < uj qo < -0.93 with x 2 = 177.1 at 68% confidence 
level. The confidence regions of Q m o and loqq are shown in figures 2 and 3. 

For the general model, the range of parameter space is Sl m0 = [0, 1], uqo = 
(-1, 0] and (3 > 0.5. The best fits to the 54 Knop sample are: < fl ma < 0.50, -1 < 
w qo < —0.23 and (3 varies in a big range with x 2 = 45.6 at 68% confidence level. 
The best fits to the 157 Riess gold sample are: < S7 m0 < 0.38, —1 < wqo < —0.57 
and (3 varies in a big range with x 2 = 177.1 at 68% confidence level. The best fits 
to the 157 Riess gold sample and WMAP data combined are: 0.21 < fi m o < 0.37, 
-1.0 < wqo < -0.75 and (3 > 0.5 with x 2 = 177.1 at 68% confidence level. From 
the above results, we see that the best fit model tends to be the A-CDM model 
with £l m0 ~ 0.3. 



5. Model-independent Results 

To construct a model independent result, we first parameterize the dark energy 
density by two parameters 32 , ^q(z) = A a + Ai(l + z) + A 2 (l + z) 2 , here Qq(z) = 
8itGpq(z)/(3Hq) and A = 1 — £l m0 — fl l0 — A\ — A 2 . The relationship between the 
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Fig. 1. The confidence contours of Q m o and ujq for the C = model. The left upper panel is 
the 68%, 95% and 99% confidence regions fitted from Knop sample. The right upper panel is the 
68%, 95%, 99% confidence regions fitted from Riess gold sample, the left lower panel is the 68%, 
95% and 99% confidence regions fitted from Riess gold and silver sample. The right lower panel 
is the 68%, 95%, 99% confidence regions fitted from Riess gold sample and WMAP data. 




Fig. 2. The 68%, 95% and 99% confidence contours of f2 m o and o;qo for f) = 1/2. The left panel 
shows fits to Knop sample and the right panel shows fits to Riess gold sample. 
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Q. _ and co_. Confidence Contour with B=0.5 
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Fig. 3. The 68%, 95% and 99% confidence contours of Q m n and ojqo for /3 = 1/2 with Riess gold 
sample and WMAP data. 




Redshitt z Redshift z 

Fig. 4. The best supernova and WMAP data fits to the polynomial model and linear model. The 
left panel shows Riess gold sample and WMAP data fits to the two parameter polynomial model, 
the light black lines are for Qq and the dark balck lines for loq, the solid lines are from the best 
fit. The right panel shows Riess gold sample and WMAP data fits to the two parameter linear 
model of uiq, the solid lines are from the best fit, the light black lines are for the linear model and 
the dark black lines are for the stable model. The dashed lines define the la boundaries. 



dark energy state of equation parameter uq and the redshift is 

_ l + z A 1 +2A 2 {l + z) 

UQ ~ 3 A + A^l + z) + A 2 (l + zf 

With the above parameteriaztion, we find that Qq f2 m and loq « — 1/3 when 
z> 1. The best fit parameters to Riess gold sample and WMAP data are f2 m0 = 
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0.40 ± 0.14, A 1 = — 3.6±| | and A 2 = 1.23±i;^| with X 2 = 174.3. By using the 
best fit parameters, we find that cjq = — 1.4 and zt = 0.37. Then we consider the 
commonly used two parameter linear model ojq(z) = ojq + wi(l + z). The Ricss 
gold sample and WMAP data give that fi m0 = 0.49±°;2o, = -2.571^8 and 
lj\ = 0.16lg^g with x 2 = 173.7. Combining the best fit parameters, it is found that 
w qo = —2.41 and zt = 0.29. Because the above linear model is divergent as z ^> 1, 
so we next consider a more stable prarmeterization of ujq(z) = ojo + u a z/ (1 + z) 33 . 
By using this parameterization, we find that the best fit parameters to Riess gold 
sample and WMAP data are O m0 = 0.47±o:i9, lu = -2.5i.4 j and w a = 3.2^2 with 
X 2 = 173.5. Therefore the turnaround redshift is zt = 0.28 and coqo = —2.5. Note 
that the cosmological constant model luq = — 1 and u a = is at the boundary of 
the la parameter space. However, the dark energy term became the dominant term 
when z ^> 1 since wo+w„ = 0.7. Therefore, we use the results from the polynomial 
parameterization only. The evolutions of Qq(z) and ujq(z) with redshift are shown 
in figure 4. 

For the model independent second order polynomial parameterization, we find 
that fi m o = 0.40 ± 0.14, wqo = —1.4 and zt = 0.37. Alam, Sahni and Starobinsky 
obtained f2 m o = 0.385 and zt = 0.39 ± 0.03 in a similar analysis 32 . Tcgmark et al. 
found that £l m0 « 0.30 ± 0.04 by using the WMAP data in combination with the 
Sloan Digital Sky Survey (SDSS) data 34 . More recently, Riess et al. showed that 
zt = 0.46 ± 0.13 from the two parameter linear model by using SNe la data only 
with the assumption that Sl m o = 0.27±0.04 4 . The above results are consistent with 
each other. For a flat universe with constant loq, we find that Sl m o = O.SOI^Qg and 
loq < —0.82. The result is consistent with our model independent results and that 
in Refs. 3, 4, 35. With those parameter values, we find that the turnaround redshift 
0.47 < zt < 0.95. For the constant pressure model /3 = 1/2, the best fits to the 
combined supernova and WMAP data are fi m o = 0.298 and ujq = —0.985 which 
result in zt = 0.65. The best parameter fits to the combined supernova and WMAP 
data for the general model analyzed in this paper are fi m o ~ 0.3, o>qo ~ —1.0 and 
(3 ~ 0.5. The turnaround redshfit is zt ~ 0.67. These results are consistent with 
the observations. In conclusion, it is shown that the general model in Ref. 20 is 
consistent with current observations and the model effectively tends to be the A- 
CDM model. Furthermore, our model independent results support the conclusion 
of dark energy metamorphosis obtained in Ref. 32. 
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